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Abstract 

In this article we prove that many Hamiltonian systems that can not be separably quantized in the 
classical approach of Robertson and Eisenhardt can be separably quantized if we extend the class of 
admissible quantizations through a suitable choice of Riemann space adapted to the Poisson geometry 
of the system. Actually, in this article we prove that for every quadratic in momenta Stackel system 
(defined on 2n dimensional Poisson manifold) for which Stackel matrix consists of monomials in 
position coordinates there exist infinitely many quantizations - parametrized by n arbitrary functions 
- that turn this system into a quantum separable Stackel system. 
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1 Introduction 

In classical mechanics the Hamiltonian equations of motion are represented by a system of nonlinear 
ode’s and are in general not integrable. A famous exception is the class of the so called Liouville 
integrable systems, i.e. those Hamiltonian systems which possess a sufficient number of global constants 
of motion in involution. In order to integrate such a system by quadratures it is necessary to find a 
distinguish orthogonal coordinates, so called separation coordinates. Once we find separation coordinates 
we can linearize equations of motion according to Hamilton-Jacobi method and then integrate them. 
Particular important class of separable systems, specially from the physical point of view, is represented 
in literature by so called Stackel systems, with Hamiltonian and all constants of motion quadratic in 
momenta. In the present paper we also restrict ourselves to such class of systems. 

This paper deals with admissible quantizations of classical Stackel systems and investigation of their 
quantum integrability and quantum separability. Surprisingly, in spite of the fact that there exists an 
extensive literature on that subject, nevertheless the foundations of the theory have been formulated in 
the early 1930’s by Robertson and Eisenhart (see the next section) and have not been changed until 
now. In their approach is considered only one particular way of quantization, which we now call natural 
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minimal quantization^ i.e. the minimal quantization generated by the metric from the kinetic part of 
the Hamiltonian of the system. One of the results of this classical theory is the so called Robertson 
condition, the fullfilment of which guarantees the quantum separability of the stationary Schrodinger 
equation generated by the corresponding quantized Hamiltonian. In consequence, according to Robertson- 
Eisenhart theory, there is only a very limited class of Stackel systems which are quantum separable. 

In this paper we broaden the theory by considering quantizations related to arbitrary metric tensor, 
not necessarily related with the Hamiltonian of the system. As a consequence of this new approach, we 
are able to formulate the following conjecture: 

For arbitrary Stackel system with all constants of motion quadratic in momenta 
there exists a family of quantizations preserving quantum separability. 

In this paper we prove that conjecture for a very large class of Stackel systems, generated by separation 
relations of the form (II3, where Stackel matrix consists of monomials in position coordinates. For any 
Stackel system from this class we construct a family of metrices for which the minimal quantization leads 
to quantum separability and commutativity of the quantized constants of motion. We want to stress, 
however, that we do not deal with spectral theory of the obtained quantum systems, as it requires a 
separate investigations. 

The paper is organized as follows. In Section[2]we briefly summarize the results of Robertson-Eisenhart 
theory of quantum separability. In Section |3] we present some fundamental facts about classical Stackel 
systems. Section|3]contains presentation of some results derived from our general theory of quantization of 
Hamiltonian systems on phase space; especially we demonstrate how to obtain the minimal quantization 
(|^ from our general theory. In Section [5] we relate quantizations of the same Hamiltonian in different 
metrics g and g (or in different Hilbert spaces Lf{Q,ojg) and L^{Q,u!g)). Essentially, this construction 
explains the origin of the quantum correction terms in the classical Hamiltonians introduced in [ij and in 
0. Section[6]is devoted to the issue of separable quantizations of Stackel systems. We construct a family 
of metric tensors which fulfill the so called generalized Robertson condition introduced in our previous 
paper Q . Using this condition we prove (Theorem [S]) that there exists an infinite family, parametrized 
by n arbitrary functions of one variable, of separable quantizations of a given Stackel system from our 
considered class. Finally, in Section [7] we address the issue of quantum integrability of Stackel systems. 
This section generalizes in an essential way the results from Q . We present the construction of commuting 
self-adjoint operators in arbitrary Hilbert spaces Lf{Q,u}g), once we have a quantum separable Stackel 
system. It also contains two illustrative examples. An invariant form of Theorem|6]is proved in Appendix. 


2 Preliminaries - legacy of Robertson and Eisenhardt 


This paper addresses the issue of separable and integrable quantizations of commuting sets of quadratic 
in momenta Hamiltonians of the form 

H {x, p) = i {x)piPg +V{x) (I) 

(throughout the whole article we apply - unless explicitly stated otherwise - the Einstein summation 
convention) defined on a cotangent bundle to some n-dimensional Riemannian manifold Q equipped with 
metric tensor g. The variables x = {xi,... ,Xn) are coordinates on Q and pi conjugate momenta (fiber 
coordinates in T*Q) while A’'^ (x) are components of a symmetric (2,0)-tensor A on Q. Note that we do 
not assume here any relation between the tensor A and the metric tensor g. The real function V{x) is 
called the potential of the Hamiltonian O- Two important partial differential equations can be associated 
with the Hamiltonian ([TJ: the Hamilton-Jacobi equation 


H [xi,... ,Xn, 


dW 
dx-i ’ 


dW 

^n. 


= a 


( 2 ) 


for the generating function W(x, a) for a canonical transformation linearizing the flow of Hamilton equa¬ 
tions 


Xi.t = 


dH 
dpi ’ 


Pi,t = 


dH 

dxF 


i = 1, 
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(here and in what follows the comma denotes the differentiation with respect to a variable) associated 
with m, and the stationary Schrodinger equation 

H'l'ix) = E^{x) (3) 

where 

H = V, + V{x) (4) 

is the Hamilton operator (quantum Hamiltonian) acting on the Hilbert space (^Q,\det dx'^ of 

square integrable (in the measure ujg = |det da;) complex functions on Q. The operators Vi are 
operators of Levi-Civita connection associated with the metric g and h is the Planck constant. One says 
then that the Hamilton operator (|T|) is the quantization of the Hamiltonian ([1} in the metric g. Note 
that the above quantization procedure is so far defined ad hoc, arbitrarily. 

An important issue related with equations ([2)) and m is the problem of their separability. We say 
that the Hamilton-Jacobi equation (I2|) is additively separable if it admits a solution 

n 

W{x,a) ='^Wt{xi,a) (5) 

i=l 

depending in a suitable manner on n additional parameters a = (oi,..., a„) (the solution ((S)) is often 
called a complete integral of ©)■ Similarly, we say that the Schrodinger equation ([3|) is multiplicatively 
separable if it admits a solution 

n 

df^x.a) = Y\_i’i{xi,a) (6) 

depending in a suitable way on 2n additional parameters a = (oi,..., a 2 n)- P- Stackel showed in Q the 
necessary and sufficient conditions for separability of © in orthogonal (with respect to A) coordinates 
(meaning that A has to be diagonal in the variables x). Assume thus that A plays the role of the 
contravariant metric (i.e. that A = G, where G = g~^) and that the metric G is diagonal in coordinates 
X. Robertson @ proved that in this case if the Hamilton-Jacobi equation (I2|) separates in the variables 
X then the Schrodinger equation ([3]) also separates provided that an additional condition, called today 
Robertson condition, is satisfied. Eisenhart in Q proved that Robertson condition is satisfied if and only 
if the Ricci tensor Rij of the metric g is diagonal. We stress again that in these works A = G and in this 
particular case the Hamilton operator Q takes the form 

JP =+ H(x) 

Robertson actually claimed in his theorem that the separability of Schrodinger equation also implies sep¬ 
arability of Hamilton-Jacobi equation: this statement is not correct if we use the definition of separability 
used by Robertson. Benenti et al in [8| completed the works of Robertson and Eisenhart by introducing 
an appropriate definition of separability of Schrodinger equation, involving 2n parameters as in 
(Robertson had no parameters in his definition of separability, a drawback not observed by Eisenhart). 
Assuming the definition of Benenti et al the theorem of Robertson becomes: 

Theorem 1 Assume that A = G and that G is diagonal in the variables Xi. The Schrodinger equation 
0) admits a separable solution 0) if and only if the Hamilton-Jacobi equation 0) admits a separable 
solution 0 and moreover if the Robertson condition 

Rij = 0 for all i j (7) 


is satisfied. 

One can show that in orthogonal coordinates 
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where are metrically contracted Christoffel symbols of g defined by 

= i = (9) 

Thus, in orthogonal coordinates the Robertson condition becomes 

diVj = 0 for j 7 ^ i (10) 

In papers @ and [ 3 ] the authors considered a quantization procedure for only one Hamiltonian and 
assumed that the underlying metric of the configuration space is defined by the tensor A in the Hamilto¬ 
nian, i.e. they assumed that A = G. Suppose now that we have n {n = dim Q) Poisson-commuting (so 
they constitute an integrable system in the sense of Liouville) Hamiltonians each of the form ([T} : 

Hr = ]^Al^PiPj+Vr{x), r = l,...,n. (11) 

A natural question one can pose is whether the corresponding quantum Hamiltonians Hr (acting in the 
Hilbert space L'^{Q,ujg) defined by the metric G = Ai) will constitute a quantum integrable systems 
i.e. whether they will commute. In the authors proved that this happens if and only if the so called 
pre-Robertson condition 

diRij - TiRij =0, j (12) 

is satisfied. Due to ([5]), this condition in orthogonal coordinates reads 

dfTg - = o, j. (13) 

Remark 2 The Robertson condition m or m implies the pre-Robertson condition m or m so 
quantum separability implies the quantum integrability, as it is in the classical case. 

The above theory describes the quantization of a Hamiltonian, or a set of Hamiltonians, of the form m 
in the case when one of the tensors Ar plays the role of the metric. However, Hamiltonians are functions 
on a phase space with no obvious metric given. In this paper we will therefore develop the theory of 
quantization of Hamiltonians of type m in Hilbert spaces L'^{Q,ujg) defined by the metric not related 
to these Hamiltonians. Let us thus pose the following question: given a separable Hamiltonian system 
consisting of n Hamiltonians of the form (El), how to find metric tensor(s) in which an appropriate 
quantization procedure turns this system into a separable and integrable quantum system? We will 
answer this question in the spirit of papers Kl where we have developed a general theory of quantizing 
Hamiltonian systems directly on the phase space; the quantization in this approach is given by an 
appropriate deformation of Poisson algebra of classical observables (real functions) on the phase space M 
to a quantum algebra. Various deformations of this algebra are related to each other by an automorphism 
S. However, in order to make this article as compact as possible, we will almost completely omit this 
general setting but use its results in the position representation, that is, we will work directly in Hilbert 
spaces L'^{Q,ujg) of the functions defined on the base manifold Q. 

3 Classical Stackel systems in separation coordinates and adapted 
Riemannian geometry 

Consider a 2n-dimensional connected Poisson manifold {M,V), where V is a non-degenerated Poisson 
tensor. An integrable system is a set of n real valued functions Hi on M in involution with respect to a 
Poisson bracket: 

{H,, Hg} := V{dH,, dHg) = 0, i, j = 1,..., n. 

The functions Hi generate n pairwise commuting Hamiltonian equations 

u^ti='PdHi, i = u £ M. (14) 

i.e. an integrable system. Let us fix a set {x,p) = {xi,... Xn,Pi, ■ ■ -Pn) of Darboux (canonical) coordinates 
on M (so that {xi^Xj} = {pi,Pj} = 0, {xi,pj} = Sij). One of the methods of solving the equations 
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(fTH) is to find a solution W{x,a) to the system of Hamilton-Jacobi equations (I2|) corresponding to the 
Hamiltonians Hi 


H^ 


f dW 

I II, • • ■ , Xm 1 



= ai, 


i = 1,..., n. 


(15) 


The solution W{x, a) is then a generating function for a canonical transformation (x, p) i—> ( 6 , a) to a new 
set of coordinates on M (with = Hi) in which the equations (1141) attain the form 


^i,tj — ^ij ; ^ 

SO that all the flows in (03 linearize in coordinates {b,a). In most cases the system of PDE’s (1151) is a 
highly nonlinear system that is very difficult to solve. However, as we mentioned in introduction, a very 
appealing situation occurs if we can find Darboux coordinates (A,/i) = {Xi,... Xn, pi, ■ ■ ■ Pn) in which 
there exists a complete integral for all the Hamilton-Jacobi equations (fTSl) of the form 


W{X,a)=Y,W^{X„a) 


(see ®) where each function Wi depend only on one canonical coordinate and in a nontrivial way on 
all parameters a = (oi,..., a„). In such a case the systems of PDE’s (fTSl) split into n uncoupled ODE’s 
for the functions Wi, which makes it possible to solve them by quadratures. The coordinates (X,p) are 
then called separation coordinates of the system (HI. 

The most convenient way to obtain separable systems is to define them directly in separation coordi¬ 
nates. It is done with the help of the so called separation relations 12|, i.e. n algebraic relations of the 
form 

ipi(Xi,fj.i,ai,...,an) =0, i = l,...,n (16) 


each depending on one pair of canonical coordinates and on parameters a^. If there exists an open dense 
set D C M on which the relations (|16D can be solved with respect to the coefficients yielding 


Ui = Hi(X, /i), i = 1,.. .,n 


then it is easy to show that the functions Hi Poisson commute (i.e. constitute a Liouville integrable 
system as defined above) and moreover that the coordinates (A, p) are separation coordinates for the 
Hamiltonians Hi. 

One of the most important classes of separable systems are the so called Stackel systems, introduced 
by P. Stackel in @ and thoroughly studied in literature (see for example 13-1^). They are generated by 
separation relations linear in Hamiltonians Hi and quadratic in canonical momenta pi. In our paper we 
restrict ourselves to a — still very general — class of Stackel systems defined by the following separation 
relations 

iJiA^+iJ 2 Ar + ---+iJ„A7" = i/,(A.)/r2 + a.(A0, i = (17) 


where 7 ^ are natural numbers such that 71 > 72 > • ■ • > 7 n = 0 (the last choice is for our convenience 
only) have no common divisor, and where fi,(Ji are some rational functions of one argument. The 
separation relations (HZ} can be written in a matrix form as 


= U, (18) 

where H = {Hi,..., Hn)’^ and U = {\fi{Xi)p\ -I- cri(Ai),..., \fn{Xn)pi + o-n(An))'^ is a Stackel vector 
and where the matrix given by 


f Xf Af ■ • • 1 \ 
V XI- ■■■ 1 ) 


is a particular Stackel matrix with functions being monomials parametrized by the natural numbers 7 ^. 
We can now take as the set D what remains of M after removing the set of points where det S'-,, = 0 as 
well as all the poles of fi and di. Solving the relations (fTSl) on D we obtain the Stackel Hamiltonians 
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r = 1,... ,n, 


( 19 ) 


with 

Ar = diag((S;\ifi(Ai), ..., (S;\nfn(An)) ( 20 ) 

being diagonal matrices with entries that are functions of A-variables only and with the potentials of the 
form 

Id" — ^ ^ (*^7 ^ri^i (Ai) T — 1, . . . , TT. 

i 

The systems of the above class, albeit not general Stackel systems, still encompass majority of the Stackel 
systems considered in literature. 

Let us now introduce some Riemannian geometry into our considerations. The specifications below 
will be motivated by the fact that our quantization procedure will be performed in appropriate (pseu¬ 
do-) Riemannian spaces. Thus, from now on we will suppose that our manifold M is a cotangent bundle 
to some pseudo-Riemannian manifold i.e. M = T*Q with Q equipped with some metric tensor g. We 
will also make three additional assumptions: 

1. The manifold {Q,g) and the Poisson structure are adapted to each other in the sense that the first 
n Darboux coordinates Xi are coordinates on Q while the remaining Darboux coordinates g,i are 
fiber coordinates. 

2. Coordinates Xi are orthogonal coordinates for the metric g i.e. g and G = g~^ are diagonal (but 
not necessarily flat) in A^. 

3. The base manifold Q is almost covered by a single, open and dense in M, chart with coordinates 

(Ai,..., Ayj). 

The matrices Ar in (fT^ can now be interpreted as (2,0)-tensors on Q that can be written as 

Ar = TrG, r = 1 ,..., n 

where are (1, l)-tensors on Q. Further, in a very special case when G = Ai the tensors are Killing 


tensors for the metric G. We will denote them as Kr, so that 

Ar = KrAi, r = l,...,n 

A particular subclass of Stackel systems (El) is then given by choosing 7 ^ = n — i. Such systems 
are called Stackel system of Benenti type (or simply Benenti systems) and are thus generated by the 
separation relations of the form 

iLiA" ^-|-i? 2 A" "^ + ■ ■ ■ + Hn =-fi{Xi)pi + o'i{Xi), i=l,...,n ( 21 ) 

It can be shown that in the Benenti case the metric tensor G = Ai has the form 

/l.=dkg(^Ll).,_i^), A.=n(A.-A,) (22) 

^ "A 

while the Killing tensors Kr are of the form 

with Pi = Pi (A) being signed symmetric polynomials (Viete polynomials) in the variables Ai,..., A„: 

p,(A) = (-!)* ^ A,, ...A,., z = l,...,n (24) 

l<si <S2<...<Si<n 
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Let us now go back to an arbitrary Stackel system of the form (HZI defined by the choice of the 
constants 71 > 72 > • • • > 7 n = 0 and the choice of functions /i, Ui. Then the tensors Ar for this system 
can be written as [13 


Ar = -XtGbj, r = l,...,n 


(25) 


where Gb.s is the corresponding Benenti metric given by (12^ 


G's,/ 



/n(A„) \ 

) 


(26) 


where Xr are some polynomial functions of the Killing tensors in (1231) and where 


(p = det 


^ Pn\—1 


V Prik-l 


Pni —k 


Puk—k 


(27) 


(where we adapt the notation pQ = 1 and pi = 0 for i < 0 or i > n) while the constants rii are those for 
which the corresponding monomials are missing in the left hand side of (flTl) (they are “holes” in 

the sequence {71 = n + fc — 1 , 72 ,-- - , 7 ra = 0} numbered from the left; k is determined from the equation 
7 i = n + k — 1). Note that if such “holes” are absent (as in Benenti case) then (p = 1. For example, if 
the left hand side of the Stackel system is HiX'^ + H 2 X + H 3 , then n = 3, fc = 2, ni = 2, 712 = 3 and the 
function (1^ becomes: 


p — det ( j — + A 1 A 2 + A 1 A 3 + A 2 + A 2 A 3 + A 3 

\ P2 Pi J 

We finish this chapter with an important remark. 

Remark 3 If fi = f and if f is a polynomial of order < n then the metric Gbj in i26\) is flat. 


4 Admissible quantizations of quadratic in momenta Hamilto¬ 
nians on pseudo-Riemannian spaces 

A usual way of quantization of a given Hamiltonian system living on a phase space M = is by 
replacing the observables of the system (i.e. real functions on the phase space of the system, written 
as functions of positions Xi and momenta pi) by self-adjoint operators acting on the Hilbert space H = 
L^(R") of square integrable complex functions on R". This is done by replacing Xi and pi in the 
observables by the non-commuting operators Xj = xj and pj = —ihd/dxj acting on L^(R"). In this 
procedure we have to agree on a certain order of non-commuting operators Xj and pj in the obtained 
operator. One usually applies the Weyl ordering that guarantees that the obtained operators will be 
self-adjoint. 

Suppose now that we want to quantize in a coordinate-free way a Hamiltonian system given on a 
phase space M = T*Q that is the cotanrant bundle to a pseudo-Riemannian manifold equipped with 
a metric tensor g. In a series of papers |9l4lll| we have developed a consistent theory of quantizing a 


Hamiltonian system directly on the phase space M through a very general procedure of quantization. 
Here we briefly sketch some parts of this construction that are important for our further considerations; 
we perform the construction in the so called position representation. 

Let us thus choose a canonical (Darboux) coordinate system {xj,pj) on M satisfying assumptions 1-3 
from the previous section. Thus, xj^ are some coordinates on Q and pj are the corresponding conjugate 
momenta. Let us also (following [l^, [l^, and (l^) introduce the operators 


Xi = Xi 


d 


1 , 


Pi = -*^ ( — + 


dxj 


(28) 
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acting on the Hilbert space H = L‘^{Q,ujg) of functions on the base manifold Q (conhguration space) 
where ojg = |det dx is a volume form defined by the metric g and where are contracted Christoffel 
symbols of the metric g. The operators (1281) are self-adjoint in % and moreover are canonical quantum 
operators as = ihSjk- 

Now, a given observable H = H(x,p) can be quantized in many different ways by applying different 
orderings to the operators x,p in H(x, p). This can be systematically done using a two-parameter family 
of automorphisms S, introduced in |ll| . acting on the space of functions on M. Any automorphism 
S from this family relates a given quantization with a Moyal quantization corresponding to our chosen 
Darboux coordinates {x,p). 

Our two-parameter family of automorphisms S is up to h^-terms given by 
S = \ + S 2 h^ + o{fi^) 

= 1 + ^ -|- aRjk)dpjdpf. + 3'Tji,dxidp^dpi^ + {2T^iT— ^jk,i)Pidpjdpf.dp^ (29) 

- Zbdp^{dx^ +T]ip,dp,)dp^{dx^ +Tl^prdpJ\+ 

(a and b are real parameters and r*^,; = 9a,,F'^.) with the inverse given formally by 

s-^ = 1 - S2h^ + o{h‘^) (30) 

Remark 4 The terms o{h^) in are at least of the fourth order in dp^ so the formulas 

are enough to calculate the action of S respectively S~^ on Hamiltonians that are up to third order in 

momenta. 


We can now introduce the following quantization procedure of a given observable H{x,p): 

1. Deformation of H(x,p) to a new function H'(x,p) = S~^H by an automorphism S from our family 

(EZl) 

2. Replacing xj and pj in H'{x,p) by the operators (E5)) . which yields the operator H'{x,p) 

3. Weyl ordering of the obtained operator. 

In short, the S'-quantization of H[x,p) in the metric g is the operator 

H = {S~^H)w{xtP) (31) 

(where W denotes the Weyl ordering) with operators x,p given by (E51) and with a chosen automorphism 
S from our two-parameter family. It can be shown that this procedure applied to any classical (real) 
observable on M yields a self-adjoint operator onH = L'^{Q,ojg). 


Remark 5 The presented procedure is invariant under the canonical change of coordinates in the sense 
that if we start from another canonical set of coordinates satisfying assumptions 1-3 from the previous 
section we obtain the quantum operator that is unitarily equivalent to H. 


Applying the above quantization procedure with the automorphism S as in (l67l) to a quadratic in 
momenta Hamiltonian 

-p'^Ap -b V(x) 


H 


yields the two-parameter family of operators (quantum Hamiltonians) on H ll| : 


(32) 


H = -y + i(l - b)A%^ - i(l - a)A^^R,j^ + V{x) (33) 

= + ?lVquant(x) + V{X) 



where is the operator of the covariant derivative of the Levi-Civita connection defined by g, Rij is the 
Ricci tensor of g and where the semicolon ; denotes the covariant derivative. The term fquant(a^) can be 
considered as a “quantum correction” to the potential V that comes from the quantization process. All 
considered in literature quantizations of quadratic in momenta Hamiltonians can be obtained by choosing 
appropriate values of a and b in (1^ . In the special case when A = G the formula (|3^ reduces to 

H = -y - i(l - a)R^ + Vix) (34) 

where R is the Ricci scalar. In the flat case (so that Rij = 0) and with & = 0 we obtain the Weyl 
quantization written in a covariant form, and (1331) and its specification (1341) attain the form 

H=-^ V, + + v{x) 

and ^ 

H = -!!-G^iV,Vj+V{x) 

respectively. 

As we see, in the general quantization scheme there appear the quantum correction term VJjuant (a;) to 
the potential V. This quantum potential is in general non-separable [l8| . so from the point of quantum 
separability the optimal choice of quantization is given by a = 6 = 1, which yields 

H = + V{x) (35) 

This quantizations is called a minimal quantization induced by the metric tensor g and Iquant(a;) = 0 in 
that case. It is exactly the a priori quantization considered by Eisenhardt, Robertson, Benenti and many 
others and described in the preliminary part above. Our theory clearly explains its origin and shows that 
this is but one of infinitely many possibilities of quantizing the Hamiltonian (HD. 

5 Minimal quantization in different metric spaces 

Our goal now is to relate two minimal quantizations induced by different metric tensors. We will need 
this in order to be able to write systems of commuting operators in various Hilbert spaces with measures 
induced by different metrics. 

Consider thus two different metric tensors g and g. As usual, we will denote their contravariant forms 
by G and G, respectively. Each of these metrics induces a minimal quantization (described in Section!?]) 
by morphisms S and S, respectively, where (cf. (l67l) with a = 6 = 1) 

*5' = 1 + [3(r;jTj^ + Rjk)dpjdpk + jk^Xidpjdpk + jk ~ ^ik,l)Pi^Pj^Pk^Pl (36) 

-idp^ (5:,. + T)ip,dp^ )dp^ (5^^ + Tl^Prdp^) + o{h'^)\ , 

and where S is given by an analogous expression with T®^ replaced by Christoffel symbols f®j, of the 
Levi-Civita connection induced by g. For a (classical) observable of the form 

H {x, p) = i {x)piPj + V (x) (37) 

by dSID, its minimal quantization with respect to g is given by 

R = (S-^H)w(x,p) = -y V,A®^V, + V(x) (38) 

and acts in L'^(Q,ujg), while its quantization with respect to g is given by a similar expression in¬ 
volving Vi (that is the covariant differentiation with respect to g) and the operators Xj = Xj and 
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Pj = These are in general two different operators, acting in two different Hilbert 

spaces: and L?{Q,ujg), respectively. The Hilbert spaces L^{Q,ijjg) and L'^{Q,ujg) are however 

isometric, with the isometry L?{Q,ujg) —>• L^{Q,ujg) given by 


jdetgl^'^^ 

= C/^' = 


|detg| 


1/4 


(39) 


where £ L'^{Q,u!g) and d' £ L‘^{Q,ujg). The isometry (15^1) induces a similarity map between operators 
in both spaces: it maps an operator F acting in L'^{Q,ujg) to the operator 


F = UFU 


-1 


acting in LF‘{Q^ojg). 


(40) 


Theorem 6 Suppose that the operator H in the Hilbert space L'^{Q,ujg) is given by 
operator UHU~^, acting in the Hilbert space L^{Q,u!g), has the form 

- hf - 

UHU-^ = -y + V{x) + h^W{x) 

with W(x) given by 

w{x) = i + 2 {r% - t%))^ 

where the subscript denotes differentiation with respect to Xi. 


Then the 


(41) 


(42) 


We will call the term W (x) the guantum correction term as it describes what happens to the operator 
transformed from L^{Q,ujg) to L^{Q,ijjg). 

Proof. One can prove this theorem by direct calculations of UHU~^. Of course 


UHU-^ = U Vj + Vix) U-^ = - — UViA^^VjU-^ + V{x) 


By using the fact 


-ffj 

dx^ 2 ^ 


after some calculations we arrive at (|4T]) - (I42I) . Alternatively, the similarity map (|40)) can be calculated 
using the automorphism SS~^. From our general theory [^- 11| it follows that quantizing the observable 
H with respect to g yields an operator that is mapped through (|40|) on the operator that we obtain by 
quantizing the observable H' = SS~^H with respect to g. This yields, that the operator (1551) attains in 
the space in L‘^{Q,u!g) the form 


UHU-^ = iS-^H')wix,p) = {S-^SS-^H)wix,p) = {S-^H)w{x,p) 


(43) 


Let us thus explicitly calculate the operator on the right hand side of (1431) . Due to (l36l) and using the 
fact that H is second order in momenta (so that the only terms in S~^ that act on H are or order up to 
h? ^ see Remark 21), after some calculations we obtain 


S-^H = H + -h^ 


1 


^ H + h^Wix) 


with 


W(x) = - 


* (r-fe - ,) + ) 

O 5* \ O \ / A \ jf- ‘'1^ Jt/ 


( 44 ) 


coinciding with W(x) in 

In Appendix we show that (1421) can be written in a covariant form as 
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(45) 




+ \A^^G^^gkr-,G^^gsi-j 


where the covariant derivatives are taken with respect to the connection V^. In what follows we will also 
need a specification of this correction term to the following situation: suppose that G = ^Gsfi (where 
u = u{x)) where the metric Gs.e is flat and suppose that G = Gsfi- Then the correction term (1^ 
attains the form 


W{x) = lU^^) 

8 V u / A 62 


(46) 


6 Separable minimal quantizations of Stackel systems 

Suppose we have a Stackel system written in arbitrary Darboux coordinates (x,p): 

Hr = ArP + Vr{x), r=l,...,n (47) 

Given a metric g we can now perform the minimal quantization of our Stackel system JUl) as described 
in the previous section. As a result we obtain n quantum Hamiltonians 

Hr = {TrGY^ Vj + Vr{x), f = 1,..., n (48) 

acting in the Hilbert space L‘^{Q,u!g), ujg = [det^l^^^ dx, where Ar = TrG. Let us rewrite the operators 
(l48t in some separation coordinates (A, p) for the classical Stackel system (ITZl) . We will always assume the 
conditions 1-3 from Section 2. This also means that g and thus G are diagonal in separation coordinates. 
Thus, since Ar are diagonal in separation coordinates, so are T^. Calculating covariant derivatives we 
obtain 


Hr = + w ( a ) 

= (^df + d}j + Vr{X) (49) 

where = {Tr)\ (no summatio^ and where T^ are metrically contracted Christoffel symbols (O. In 
orthogonal coordinates they read y 


1 di det G 

* “ 2 detG ^ 

The next theorem, proved in [^, follows directly from (l49l) . 

Theorem 7 The necessary and sufficient condition for quantum separability of operators Hr takes the 
form 

Ei =Ei{Xi) or djEi = 0, j i (50) 

where 

„ _ 

We will call the condition (1501) the generalized Robertson condition. Indeed, due to (1001) , the operators 
(HOt can then be written as 

Hr = -\h^{S-XMX^){df + E,iXYd.) + {SfXa.iXY, r = l,...,n (51) 
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and then application of the Stackel matrix to the system of eigenvalue problems for (IFTl) 


5^ 






= 




\En'^/ 


( 52 ) 


separates (1521) to n one-dimensional eigenvalue problems 


{E,X'J^+E2XT+--- + En)M^,) = --h‘^MX^) 


(Pi>t{Xi) 




dipiiXi) 


dX^i ' dA, 

called separation equations or quantum separable relations, so that 

n 

^'(Ai, . . . ,Xn,C,E) = W_'lpi{Xi,C2i-lC2i,E) 


+ cF^{Xi)'ipi{Xi), i = 

(53) 


is a common, multiplicatively separable solution of stationary Schrddinger equations for all satisfying 
the definition of separability from . The constants Ei are unspecified unless some boundary conditions 
are imposed while C 2 i-iC 2 i are integration constants originating during the process of solving equation 
i in (1551) : there are 2n of them in total. In the case G = Ai (or, in general, G equal to any A^) Tr 
are Killing tensors of g so in A-coordinates diTr^^ = 0. In consequence the condition (150)) reduces to the 
Robertson condition for quantum separability 0 or (nni). 

In [l^ we proved that for the case G = Ai the only class of Stackel systems m for which the 
Robertson condition dm is satisfied is the Benenti class where 


1 /'(A.) 

2 MX,) 


(54) 


For all other choices of % in (ITTl) this condition fails. In Q we investigated the more general case when G 
is not one of the tensors Ar in dm but is a flat metric from the Benenti class (l26l) . We showed that also 
in this case the only class of Stackel systems dm that is quantum separable is again the Benenti class. 
It means that in order to achieve quantum separability of an arbitrary Stackel system of the type dm 
we have to consider a broader class of admissible metric tensors g used in the quantization procedure. 

Consider thus a Stackel system dm defined by some fixed choice of 71 > 72 > • • • > 7 n = 0 and the 
choice of We will now search for the metric G that satisfies the generalized Robertson condition 

(|50)) for this Stackel system. Due to the structure (|25|) of Ar we look for G in the form 


G = u-\X)GB,e 


(55) 


where Gb.b is the Benenti metric given by (1261) with n arbitrary functions 0i{Xi) and where u is some 
function on Q. Albeit this choice is by no means the most general one it will prove to be sufficiently 
general. The tensors Tr become in this case 


Tr = —XrGB,fgB,e 


where g} is again given by (1571) and where as usual gB,e = Gg^g. Plugging this into (1551) we get 


d,TM 


MM 

B'i (Aj ) 


= I,...,n 


(56) 


where Ki are arbitrary functions of one variable (the right hand side is just a convenient for us way of 
writing an arbitrary function of Ai). Since for (1551) 


Pi 


(PB,e)i + 



d,u 

u 
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with (r^.e), being the metrically contracted Christoffel symbols for the metric GB,e 
takes the form 

nd,u djip _ k[{K) //(A,) 

2 u ip Ki{Xi) 2 9i{Xi) fi{Xi)' 

which has a solution 


u = 


n 

2=1 


m 


J i 


In order to receive a solution as simple as possible we choose Ki so that 


, the formula (155)) 


(57) 



= 1 


(notice that 6 i are still arbitrary) yielding (1571) in the form u = (p^ . Thus, we have proved 

Theorem 8 Suppose 0i,i = 1, ..., n are n arbitrary functions of one variable. Then applying the proce¬ 
dure of minimal quantization, with the metric tensor 

g = T" 9 B,e (58) 


where gsp = with Gb ,0 given by 


Gb.o = diag 


V 


0n(A„)\ 

A„ J 


(59) 


to the Stdckel system we obtain a quantum separable system with the separation equations of 
the form 


iEiXT+E2X]^ + ---+E„)MX,) 


where i = 1,... ,n. 




df'ifi{Xi) 

dXl 



1 

2 e,{x,)) 


difijXj) 

dXi 


-\- (r iifXi^tf i{Xi), 

(60) 


The metric g in (1581) is a conformal deformation of the Benenti metric gB,e- Thus, there exists an 
infinite family of separable quantizations of a Stackel system dni) parametrized by n arbitrary functions 
9i of one variable: any Stackel system (HZ} can be separably quantized in the conformally deformed 
metric (1551) (note that this metric is conformally flat in the case when gB,e is flat). Moreover, since for 
the Benenti class (p = 1, any Stackel system from the Benenti class (1511) can be separably quantized in 
any metric of Benenti class (1591) . including the subclass of flat metrics. 


7 Quantum integrability of Stackel systems in arbitrary Hilbert 
spaces 

We remind the reader that in the authors derived the necessary and sufficient condition for com¬ 
mutativity of quantum Hamiltonians of the form dill) (and with Ai = G) called the pre-Robertson 
condition (EH) or m, which took the form 

dfTj - = o, ^ ^ j. ( 6 i) 

In our case, when G is not related with any Ar, analogous calculations lead to the following necessary 
and sufficient condition for commutativity of Hr which we call the generalized pre-Robertson condition 

0 : 

=0, i^ j. (62) 

Assume that we have a Stackel system Hr, r = 1,... ,n oi the form dn}. Let us perform the procedure 
of minimal quantization of this system in the metric G given by d58|) , as described in the previous section. 
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We obtain then the quantum separable system consisting of n operators acting on the Hilbert space 
L‘^{Q,ujg), ujg = \Aetg\^^^ dX. Since the generalized Robertson condition (15(11) implies the generalized 
pre-Robertson condition (15^ we conclude that this system is also quantum integrable: Hr,Hs = 0. 


Using Theorem|6]we are able to write operators Hr in another metric g i.e. in the Hilbert space L^{Q, ujg) 
which yields new quantum operators Hr-: r = 1,..., n that constitute again quantum integrable (but not 
necessarily quantum separable) system. Due to the theory developed in Section [5] we know, that we can 
equally well take the classical Hamiltonians Hr amended by quantum correction terms, i.e. the functions 
Hr + h^Wr with Wr given by (or equivalently by (05])) and minimally quantize them in the metric 


g as this will yield the same quantum integrable system Hr-, r = 1 ,..., n. 

In [l^ we demonstrated that any Stackel system of the class (13 can be constructed by an appropriate 
Stackel transform of a suitably chosen flat Stackel system from Benenti class. Moreover, in [20[ we 
explicitly constructed flat coordinates for any flat Stackel system. Therefore we are able to write down 
our original Stackel system Hr, r = 1, ... ,n in flat coordinates of the metric g of the form (l26l) { g is 
flat as soon the conditions in Remark O are satisfied). In this specific case, if we apply the standard 
Weyl quantization to the Stackel system Hr + h^Wr (i.e. our original system amended by the quantum 
correction terms h^Wr) we will obtain a quantum integrable system. One can also say, alternatively, that 
if we want to avoid quantum correction terms, we should quantize the original system Hr, r = 1 ,... ,n 
not by Weyl quantization but by minimal quantization in a suitably chosen conformally flat metric G. 

In papers [ij and Q the authors presented some ad hoc calculations generating quantum correction 
terms that guarantee integrability of quantum systems obtained through Weyl quantization of some 
Hamiltonian systems. Our theory shows how to construct these quantum correction terms in a systematic 
way (albeit within the class of Stackel systems, not considered in i-i). We will illustrate this on two 
examples below. It is important to stress that the presented systems cannot be separably quantized in 
the frame of the classical Robertson-Eisenhart formalism. 


Example 9 Consider the Stackel system for n = 3 given by the separation relations of the form: 

HlX^,+H2X, + H3 = ^X^g^,+^t, i = 1,2,3 (63) 

so that 7 i = 3, 72 = 1 and 73 = 0 and with fi{Xi) = Xi and cri(Ai) = A^. In this case ip = pi(A) = 
— (Ai -|- A 2 -l- A 3 ). Consider also the corresponding metric Gbj given by 126]) . This metrie is fiat, by 
Remark\^ In the coordinates xi,X 2 ,X 3 defined through (cf. \2^) 


the metric Gs.f reads 


pi = — (Ai -I- A 2 -I- A 3 ) = xi 

1 2 

P2 = A 1 A 2 -|- A 1 A 3 -|- A 2 A 3 = a ;2 -l- — 

P3 = —A 1 A 2 A 3 = 

/ 0 1 0 
Gs.f =10 0 
\ 0 0 1 


(64) 


(65) 


so (p = Xi in Xi-coordinates and Xi are flat non-orthogonal coordinates for Gbj- Solving the relations 
in with respect to the Hamiltonians Hi and passing to the variables Xi we receive Hr = A^fpiPj -\-Vr{x) 
where pi are momenta conjugate to Xi and where the tensors Ar have the form 


j 


_ J_ 

Xi 


i j 

f 1 

ixi - 

4 x-]_ 


Ali = 

, -ST 

0 


, ^2 = 


-X2 

-2^3 

1 

^ 0 

0 


' 1 

^ 0 

-\X3 

!^i - ft/ 


AI 3 = 


1 ^ 2 . 
4 XI 

-ix3 


1 

4 xi 
4*^3 


- 5 X 3 

-7X1X3 


-iXiX3 M+X2 + i^J 
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with the corresponding rational potentials 


= -i« +5 


Vsix) 


16 


xixl 


1 X 2 ^ 
4 xi 


From our theory it follows that we can perform a separable quantization of this system in the conformally 
flat metric G = ^Gbj (which means that we choose 6i = fl) with u = We obtain three 

commuting operators 

Hr = + Vr{x) (66) 

(where is the connection defined by G),acting in the Hilbert space L^{Q,ujg) = L^{Q,\xi\dx) (ujg = 
|det 51^^^ dx = dx = \xi\dx). In the separation coordinates (A,/i) the separation equations S6(M for 

Hr attain the form 


{Ei\\ + E2\i + 


lfe 2 A 


1 dil}i{\i) \ 

2 d\, ) 


i = l,2,3 


(67) 


Let us now rewrite our operators (EW in the Hilbert space L^{Q,ujfl) = L^{Q,dx) (uig = Idet^l^^^ dx = dx) 
with the flat metric G = Gbj- From our theory it follows that a suitable way to do it is to quantize our 
Hamiltonians Hr directly in the metric G after amending them by the quantum correction terms Wflx) 
given by 

Wi = 0, W 2 = ^3 = --Z — 

8x1 8 a;i 

One can check by direct calculations that the operators 


Hr = --hfdiAfidj + h'^Wr{x) + Vr{x), v = 1, ..., n 


( 68 ) 


(the coordinates Xi are flat for g = Gbj so = di = d/dxi) do indeed commute, thus constituting 
a quantum integrable system. The operators lEW are however not quantum separable, contrary to the 
operators (EE), but are R-separable. It means that in separation coordinates 


Hr'^iX) = Er'i>{\), ^(A) = C/(A)4'(A) = (A^ + A^ + 


and ip{X'') solves m- 

Example 10 In our second example we consider the following Stdckel system 

HiXl+H2X^,+H3 = ^X^fj.( + Xt, i= 1,2,3 (69) 


so that this time 71 = 3, 72 = 2 and 73 = 0 but still with fi{Xi) = Xi and cri{Xi) = Xf. In this case 
ip = P 2 (A) = A 1 A 2 + A 1 A 3 + A 2 A 3 . We consider again the same metric Gbj with the same flat coordinates 
Xi given by This time the tensors Ar have the form 



^ —1 —\x\ 0 \ 


/ -Xi -Axi^+X2 0 \ 

P2(x) 

-Axi X2 5X3 

’ "^2 - /X 

P2(X) 

-Axi^+X2 X2X1 5X1X3 


y 0 Ax 3 -Xi y 


\ 0 5X1X3 -fxi^ +X 2 J 


/ X-i^ 

1 2 

2 ^3 

-A (xi^ + 4x2) X3 \ 

. _ 1 
^ 4 p 2 (a;) 

1 2 

2 ^3 

jXif^x F 

-A X3 (xi^ +4X2X1 +2X3^) 


\ (a;i^ + 40:2) 0:3 X3 (xi^ + 4 x 2 X 1 + 2x3^) A xi"^ + 2xi^X2 + 4 x 2 ^ + X3^xi J 
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where P 2 {x) = X 2 + \x\, while the potentials are 


V 2 {x) = - ^ ^1^4 , 0.2 


4 ^ 2 ( 2 ;) V4 


Xj + 2xfx2 + XiX^ + 4^2 


Vsix) = 


16 p 2 (x) 


This time we perform a separable quantization in the conformally flat metric G = ^Gbj with u = qp’!'^ = 
( 2^2 + obtain again three commuting operators 


Hr = --h^V^A'-fVj+Vrix) 


(70) 


(where is the connection defined by G),acting in the Hilbert space L'^{Q, |p 2 (a;)| dx), while the separation 
equations \6(A) for Hr become 


{EiX'( + E 2 PI + E3)ifi{Xi) — --h^ ( Xt — 


1 fc2 /^x d^ifflXi) ^ 1 dipi{Xi) 

2 


2 dXi 


+ XtfjflX,), I = 1,2,3 (71) 


with the same right hand side as in the previous example. Rewriting our operators in the Hilbert 
space L^{Q,dx) with quantization defined by the flat metric G = Gbj leads to the following correction 
terms Wi{x) 


Wi = 


1 


(5x^ - 4 x 2 ) 


Again, the operators 


I 6 pl{x) 

= 3 ^ 7-f - 20xw.) 

^3 = X (a^i + 8 x^X 2 + 13X1X3 + 16x1X2 + 4x2X3) 

128 p^(x) 

Hr = --hfdiA^Jdj + fl^Wrix) + Vr{x), r = 1,..., n 


(72) 


commute, as it can be checked for example in Maple. Operators are R-separable and in separation 
coordinates 

Hr^{X) = Er^{X), ’I'(A) = t/(A)4'(A) = (A^A^ + A^A^ + X^X^)^fj{X^)i:{X^)-ilj{X^), 

where V'(A*) solves m- 


8 Appendix 

We sketch here the proof of the fact that formulas (l42l) and (1^ are equivalent. We want to demonstrate 
that 


w(x) = - A^^ + 2 {a^= {r% - t%)) ^ 


(73) 


or, equivalently 


W(x) = 


1 


— tpfc _ pfe I 1 Aij 1 4- i /pfc pi _ pfc pi \ 


(74) 


coincides with 
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(75) 


W{x) = i + \A^^G^^gkr-,G^^gsi-: 

where the covariant differentiation is taken with respect to the metric g. To this end, denote by M) the 
quotient of gij and gij: 

(76) 

(as such, it is a (1, l)-tensor), which yields 

The Christoffel symbols T*^, can now be expressed through T*^ in the following way 


i9ij,k + gik,j - g]k,i) 


— {^i,k9sj + Migsj,k + M^jgsk + Migskj — MJ^gsk — M-gsk,i) 


yielding 


Using 


rjfc = fjfc + + ^{M-^)l.Miyg^^gsk - ^{M-^)lMl,f^g,k 

- ^-{M-^)iM^r%k,i + lr9jk,s. 


0 — 9jk-,s — 9jk,s SrifeT g gjnTi.g, 


(where ; denotes the covariant differentiation with respect to g) we receive 


= ^k + liM-^TrMJ.k + - l{M-^)lMJ,if^g,k. 


In particular 


Moreover 


= ^Jk + ^{M-^YrMl^y 




Inserting all this into (l7^ we obtain 




that due to (1751) coincides with (1751) . 
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